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ABSTRACT 

We present the results of a survey for Lya emitters at z w 5.7 based on optical narrow-band (A c = 8150 
A and A A = 120 A), and broad-band (B, Rc, Ic-, an d z') observations of the field surrounding the high 
redshift quasar, SDSSp J104433.04-012522.2, on the 8.2 m Subaru Telescope with the Subaru Prime 
Focus Camera, Suprime-Cam. This survey covers a sky area of ~ 720 arcmin 2 and a co-moving volume 
of ~ 2 x 10 5 7 Mpc 3 . We have found 20 candidates of Lya emitters at z ~ 5.7 with Az w 0.1. Two of 
them have been confirmed star-forming galaxies at z = 5.655 and z = 5.687 from our follow-up optical 
spectroscopy. We discuss star-formation properties of the 20 objects from a statistical point of view. 
Our survey leads to a new estimate of the star formation rate density at z » 5.7, ~ 1.2 x 10 _3 /io.7-^© 
yr -1 Mpc -3 . 

Subject headings: cosmology: observations — cosmology: early universe — galaxies: formation — 
galaxies: evolution 
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wide-field CCD cameras on 8 - 10 m class telescopes. The 
Suprime-Cam mounted at the prime focus of the 8.2 m 
Subaru Telescope at Mauna Kea (Miyazaki et al. 2002) 
provides a unique opportunity for wide-field (a 34' x 27' 
field of view), narrowband imaging surveys for emission- 
line objects at high redshift. The efficiency of this in- 
strument is higher by a factor of 30 than typical imagers 
mounted on the other 8 - 10m class telescopes. 5 In this 
paper, we present results of our imaging survey using this 
instrument to investigate observational properties of star- 
forming objects at z w 5.7 (see also Ajiki et al. 2002; 
Shioya et al. 2002; Taniguchi et al. 2003). Our survey 
field includes the quasar SDSSp J104433.04-012502.2 at 
z = 5.74 6 . 

We adopt a flat universe with f2 ma tter = 0.3, f^A = 0.7, 
and ho.r — 1 where ho.7 — Hq/(7Q km s _1 Mpc -1 ). 
Throughout this paper, magnitudes are given in the AB 
system. 

2. OBSERVATIONS AND DATA REDUCTION 

2.1. Observations 

We have carried out a very deep optical imaging survey 
for faint LAEs in the field surrounding the quasar SDSSp 
J104433.04-012502.2 at redshift 5.74 (Fan et al. 2000; 
Djorgovski et al. 2001; Goodrich et al. 2001), using the 
Suprime-Cam (Miyazaki et al. 2002) on the 8.2 m Subaru 
Telescope (Kaifu 1998) on Mauna Kea. The Suprime-Cam 
consists of ten 2k x 4k CCD chips and provides a very wide 
field of view, 34' x 27' (0.202 arcsec pixel" 1 ). In this sur- 
vey, we used a narrow-passband filter, NB816, centered on 

^^Based on data collected at Subaru Telescope, which is operated by the National Astronomical Observatory of Japan. 

2 Based on data collected at the W. M. Keck Observatory, which is operated as a scientific partnership among the California Institute of 
Technology, the University of California, and the National Aeronautics and Space Administration. 

3 Astronomical Institute, Graduate School of Science, Tohoku University, Aramaki, Aoba, Sendai 980-8578, Japan 

4 Subaru Telescope, National Astronomical Observatory of Japan, 650 North A'ohoku Place, Hilo, HI 96720 

5 Note that several wide-field cameras are also available on 4m class telescopes such as the Hale 5m telescope, KPNO and CTIO 4m telescopes, 
and CFHT (e.g. Roads et al. 2000; Stiavelli et al. 2001). 

6 The discovery redshift was z = 5.8 (Fan et al. 2000). Since, however, the subsequent spectroscopic observations suggested a bit lower 
redshift; z = 5.73 (Djorgovski et al. 2001) and z = 5.745 (Goodrich et al. 2001), we adopt z = 5.74 in this paper. 
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1. INTRODUCTION 

Lymana emitters (LAEs) are one of the important pop- 
ulations of star- forming galaxies in the early universe (e.g., 
Partridge & Peebles 1967; Larson 1974; Meier 1976). 
Therefore, the search for LAEs at high redshift is essen- 
tially important for the understanding of both the his- 
tory of galaxies and the origin of cosmic reionization (e.g., 
Loeb & Barkana 2001). Recent optical surveys with a nar- 
rowband filter have been finding LAEs beyond redshift 5 
thanks to the great observational capability of 8 - 10 m 
class telescopes (e.g., Hu, Cowie, & McMahon 1998; Hu et 
al. 2002; Rhoads et al. 2003; Kodaira et al. 2003). 

However, such surveys have the following two technical 
limitations. The first limitation comes from strong OH air- 
glow emission lines because they dominate at wavelengths 
longer than 7000 A (e.g., Stockton 1999). Although there 
can be seen several gaps at which there is little OH emis- 
sion, this prevents us from finding very faint galaxies at 
high redshift. The well-defined gaps appear around A w 
7110 A, 8160 A, and 9210 A. These gaps enable us to 
search for LAEs at z w 4.8 (Ouchi et al. 2003; see also 
Rhoads et al. 2000), 5.7 (Hu et al. 1998; Rhoads & Mal- 
htora 2001; Ajiki et al. 2002; Taniguchi et al. 2003), and 
6.6 (Hu ct al. 2002; Kodaira et al. 2003). The second 
limitation is that survey volumes are so small because of 
narrow filter band widths (e.g., ~ 100 A). 

Although there are the above two limitations, recent 
successful results encourage us to promote a new deep 
survey with a narrowband filter. In order to gain survey 
volumes and to reach faint limiting magnitude, we need 
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8150 A with a passband of AAfwhm = 120 A; the wave- 
length corresponds to a redshift of 5.65 - 5.75 for Lya 
emission 7 . We also used broad-passband filters, B, Rc, 
Ic, and z'. The total-response (filter, optics, atmosphere 
transmission, and CCD sensitivity are taken into account) 
curve of the each filter band used in our observations is 
shown in Figure 1. 

A summary of the imaging observations is given in Table 
1. All observations were done under photometric condi- 
tions, and the seeing was between Of.' 7 and l'/4 during the 
observing run. Photometric and spectrophotometric stan- 
dard stars used in the flux calibration are SA101 (Landolt 
1992) for the B, R c , and I c data, and CD 50, GD 108 
(Oke 1990), and PC 1034+001 (Massey et al. 1988) for 
the NB816 data. The z' data were calibrated by using the 
z' magnitude of SDSSp J104433. 04-012502. 2 (Fan et al. 
2000). Since any quasar is a potentially variable object, 
the photometric calibration of the z' data may be more 
unreliable than those of the other-band data. Then, using 
colors of stars obtained from spectral energy distributions 
(SEDs) of various spectral types of stars (Gunn & Stryker 
1983), we have made a small amount of correction (+0.15 
mag) for the z' data. The individual CCD data were re- 
duced and combined using IRAF and the mosaic-CCD 
data reduction software developed by Yagi et al. (2002). 

The total size of the reduced field is 32(08 x 25(23 « 807 
arcmin 2 . After masking the regions contaminated by 
fringes and bright stars, our actual survey area is w 720 
arcmin 2 . The final image of NB816 is shown in Fig- 
ure 2 together with actual survey area enclosed by solid 
line. The volume probed by the NB816 imaging has (co- 
moving) transverse dimensions of 4.0 xl0 3 ft.Q~ 2 Mpc 2 , and 
the FWHM of the filter correspond to a co-moving depth 
along the line of sight of 44 \ Mpc (z m j n w 5.65 and 
Zmax ~ 5.75; note that the transmission curve of our 
NB816 filter has a Gaussian- like shape). Therefore, a to- 
tal volume of 1.8 x lO 5 ^ 3 Mpc 3 is probed in our NB816 
image. 

2.2. Source Detection and Photometry 

Source detection and photometry were performed using 
SExtractor version 2.2.2 (Bcrtin & Arnouts 1996) double 
image mode. A source is selected as a 13-pixel connection 
above the 2 a noise level on the NB816 image. Photom- 
etry is performed with a 2.8 arcsec diameter aperture for 
each band image of which image sizes are matched with 
the i?c _ band image (1.4 arcsec). The limiting magnitudes 
are NB816 = 26.0, B = 26.6, R c = 26.2, I c = 25.9, and 
z' = 25.3 for a 3a detection with a 2.8 arcsec diameter 
aperture. In the above source detection, we find ~ 49, 000 
sources down to NB816 = 26.0. 

3. RESULTS 

3.1. Selection of NB816-Excess Objects 

Since the central wavelength of the Iq filter is bluer than 
that of the NB816 filter, we calculate a continuum magni- 
tude at A = 8150 A that we refer to as the "Iz continuum" , 
using a linear combination (fj z = 0.76// c +0.24/ z /) of the 
Ic and z' flux densities; a 3 a limit of Iz ~ 26.0 in a 2.8 

7 In Ajiki et al. (2002), the central wavelength of the NB816 filter i 
The redshift range covered by this filter is given as 5.663 - 5.762. Thes 
discussion in Ajiki et al (2002). 



arcsec diameter aperture. This enables us to sample more 
accurately the continuum at the same effective wavelength 
as that of the NB816 filter. 

We select NB816-excess objects using the following cri- 
teria; 

NB816 < 25.0, (1) 
Iz - NB816 > 0.9, (2) 

where Iz - NB816 > 0.9 correspond to EW ohs £ 180 A. 
In Figure 3, we show the diagram between Iz — NB816 
and NB816 for the objects in the NB816-se\ected catalog. 
There are 62 sources which satisfy the above two criteria. 

3.2. Selection of LAE Candidates 

Our main aim in the present survey is to find strong 
LAEs at z « 5.7. However, strong emission-line sources 
at lower redshift may be also found in our survey; e.g., [O 
n]A3727 sources atzw 1.2, H/3 A4861 sources atz« 0.68, 
[O lli]A5007 sources at z « 0.63, Ha+[N n] A6583 sources 
at z « 0.24 (see Fujita et al. 2003b), and so on. In order 
to select LAEs at z w 5.7 from our emission-line objects, 
we apply the following criteria to all emitters selected in 
the previous section; 

B > 26.6, (3) 
R C -I C > 1-8 for I c < 24.8, (4) 
R c > 26.6 for I c > 24.8. (5) 

The first criterion means that objects at z « 5.7 must be 
undetected at 3 a noise level in our B band image. Be- 
cause the continuum break at Ao = 912 A (Lyman Break) 
shifts to A bs ~ 6200 A, galaxies at z « 5.7 are completely 
damped in _B-band (A max « 5000 A) . The second criterion 
is led as follows. We show the diagram between Rc — Ic 
and Ic — z' for all the objects detected in the broad-band 
images in the right panel of Figure 4. We also show color 
evolutions of model galaxies with the SEDs typical to SB 
(young starburst), Irr, Scd, Sbc, and E as a function of red- 
shift (left panel). As for the generation of SEDs of model 
galaxies, we use GISSEL96 (Bruzual & Chariot 1993). The 
star formation rate (SFR) of model galaxies is proportional 
to exp(— t/r), where we use r = 1 Gyr. Changing the age 
of model galaxies, we generate the SEDs of above five type 
galaxies (t = 1, 2, 3, 4, and 8 Gyr for SB, Irr, Scd, Sbc, 
and E, respectively). We use the average optical depth 
derived by Madau et al. (1996) to estimate the absorp- 
tion by intergalactic neutral hydrogen. The model results 
show that LAEs at z <~ 5.7 may have Rc — Ic J; 1-8 and 
Ic — z' ^> 1.0 although the latter color constraint appears 
not so severe. This figure also implies that lower- z emit- 
ters could have Rc — Ic < 1-8. However, the majority of 
our emitters are too faint (Ic > 24.8) to be applied by 
this color criterion. For such faint sources, we set another 
criterion; LAE candidates should be undetected above 2 a 
in the i?c _ band image (the third criterion). By applying 
the above criteria to our 62 emitters, we select 20 LAE 
candidates at z ~ 5.7 and none of them has Ic < 24.8. 

given as A ccn t C r = 8160 A. This should be read as A con t C r = 8150 A. 
values should be read as 5.65 - 5.75. These changes do not affect any 
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Their positions and photometric properties are given in 
Table 2. It is noted that all of our LAE candidates are 
undetected above 3<r level in the B- and i?c-band images 
(i.e., B > 26.6 and R c > 26.2). In Figure 5, we show 
the B, Rq, Ici NB816, and z' images (upper panel) and 
contour maps (middle panel) of the LAE candidates at 
z « 5.7. The SED is also shown in the lower panel for 
each object. We note that the two objects, No. 17 and No. 
18, have been confirmed as galaxies at z = 5.655 ± 0.002 
(Taniguchi et al. 2003) and z = 5.687 ±0.002 (Ajiki et al. 
2002) by spectroscopic observations (see Section 4.6). 

3.3. Contamination by Low-z Emitters 

There is a possibility that some low-z emitters are in- 
cluded in our final sample of the 20 LAE candidates. Since 
this affects our later discussion, we examine this possibility 
more quantitatively. 

All the LAE candidates are very strong emission-line 
sources because they have EW b s ^ 180 A. If they were 
lower-z emitters, they would be star-forming galaxies at 
z = 0.24, 0.63, or 1.19 corresponding to Ha, [O m]A5007, 
or [O n]A3727, respectively. They could be very strong 
emission-line galaxies because our observed equivalent 
width limit corresponds to EW{Ra) ;> 145 A, EW([0 in]) 

;> 110 A, and EW([0 n]) £ 82 A, respectively. For the 
above three cases, we can nominally estimate the SFR for 
the 20 NB816-excess objects (i.e., the LAE candidates.). 
Then we can estimate their apparent B and Rc magni- 
tudes based on galaxy models generated by using GIS- 
SEL96 (Bruzual & Chariot 1993). If these magnitudes are 
bright enough to be detected in our B and Rc images, we 
can reject the possibility that they are low-z emitters. 

First, we estimate the nominal SFR for the above three 
cases; i.e., Ha, [O n], and [O ill] emitters. 

(1) The case of Ha emitters: We cannot estimate a pure 

Ha line flux for each source from our data because 
we have no information on the contribution of [N 
n]A6583 emission to the observed flux. Therefore, 
for simplicity, we adopt the average line flux ratio 
for star-forming galaxies, /(Ha)//([N n]) = 2.3, ob- 
tained by Kennicutt (1992), and then estimate the 
pure Ha flux from our data. 

We convert the Ha luminosity to the SFR for each 
object, using the following relation; 

SFR(Ra) = 7.9 x 10~ 42 £(Ha) A^yr" 1 (6) 

where, -L(Ha) is the Ha luminosity in units of ergs 
s~ 4 (Kennicutt 1998). In this estimate, the solar 
metallicity (Z = Z Q = 0.02) and the Salpeter initial 
mass function with 0.1 - 100 M Q are assumed. 

(2) The case of [O II] emitters: We use the following re- 

lation; 

SFR([0 ii]) = 8.8 x 10~ 42 L([O n]) A^yr" 1 (7) 

where L([0 n]) is the [O n] luminosity in units of 
ergs s _1 . In this estimate, we use equation (6) to- 
gether with [O n] / Ha = 0.9 (Hippeline et al. 2003). 



(3) The case of [O III] emitters: We use the following 
relation; 

SFR([0 in]) = 1.0 x 10~ 41 L([O in]) Moyr- 1 (8) 

where L([0 in]) is the [O in] luminosity in units 
of ergs s _1 . In this estimate, we use equation (7) 
together with [O m]/[0 n] = 0.9 (Hippeline et al. 
2003). 

Next, we estimate the expected apparent B and Rc 
magnitudes of model galaxies with the SFRs estimated 
above, given that objects are located at the corresponding 
rcdshifts. In order to estimate apparent B and Rc magni- 
tudes, we use GISSEL96 assuming that these galaxies have 
been forming stars constantly for the recent past, 0.01, 0.1, 
or 1 Gyr and have solar metallicity (Z = Z Q = 0.02). The 
results are shown in the three left panels of Figure 6 for 
the case of Ha, [O in], and [O n] emitters. 

In the above estimates, we do not include the effect of 
internal extinction. However, since the internal extinction 
affects both the line luminosity and the apparent magni- 
tude, we examine how the extinction affects our results. 
Since it is known that typical starburst galaxies have the 
color excesses for the nebular gas, E(B — V) gas = - 
0.6 (e.g., Calzetti et al. 2000; Cardiel et al. 2003), we 
also estimate expected apparent B and Rc magnitudes 
for E(B — V) gas = 0.3 and 0.6. In these estimates, 
we adopt the extinction law obtained by Calzetti et al. 
(2000). The results for the cases of E(B - F) gas = 0.3 and 
E(B — V) gas = 0.6 are shown in the second column and 
the right column of Figure 6, respectively. 

In Table 3, we give a summary of our results. As shown 
in Figure 6, if our LAE candidates were low-z emitters, 
they could be easily detected in our B and Rc images. 
Therefore, we can rule out the possibility that our LAE 
candidates are low-z emitters. This conclusion is hold even 
if the internal extinction works up to E{B — V) gas = 0.6. 
The reason why the effect of the extinction does not play 
an important role is that the obscuration by dust is more 
effective for the line emission from nebular gas than that 
for the stellar-continuum emission i.e., 

E(B - V) star = 0.UE(B - V) gas (9) 

(Calzetti et al. 2000). One remaining possibility seems 
that some candidates may be Ha emitters at z « 0.24 if 
they experience a short-timescale star formation event (the 
duration must be shorter than ^0.1 Gyr). However, it is 
unlikely that there arc numerous such galaxies at z = 0.24. 
Therefore, we conclude that our LAE sample is free from 
the contamination of low-z emitters. 

3.4. Star Formers vs. Active Galactic Nuclei 

Previous narrowband imaging surveys for high-z 
emission-line sources resulted in the discovery of star- 
forming galaxies (e.g., Hu et al. 1998, 2002; Kodaira 
et al. 2003; Rhoads et al. 2003). However, it is uncer- 
tain whether our LAE candidates are genuine star-forming 
galaxies or active galactic nuclei, especially quasars, al- 
though our two LAEs are star forming galaxies (Ajiki et 
al. 2002; Taniguchi et al. 2003). Our survey is biased 
to detect very strong emission-line objects (EW b s J> 180 
A). Therefore, if observed equivalent widths of Lya+N 
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v emission lines are smaller than the EW limit, quasars 
are not selected as LAE candidates. For example, SDSSp 
J104433.04-012502.2 at z = 5.74 (Fan et al. 2000) and 
SDSS J114816. 64+525150.2 at z = 6.43 (Fan ct al. 2003) 
have EW hs ~ 180 A 8 . However, typical quasars have 
EW - 70 A (e.g., Fan et al. 2001). In fact, some very 
high-z quasars around z ~ 6 have EWq ~ 70 A (e.g., Fan 
et al. 2002, 2003). If such quasars at z « 5.7 are present 
in our survey field, they could be selected as LAE can- 
didates because EW b s exceeds 450 A. Although broader 
line widths of Lya+N v emission lines make it difficult 
to detect high-z quasars in any narrowband imaging sur- 
veys, it seems important to estimate a probable fraction 
of quasars in our LAE sample. 

For this purpose, we estimate the expected number of 
quasars, iVquasar, in our survey volume. In the estimate, 
we assume that the number density of quasars at z w 5.7 
follows the maximum likelihood solution derived by Fan 
et al. (2001) based on their data of 39 luminous quasars 
at z = 3.6 -5.0; 



log$(z,< M1450) 



(-7.31 ±0.19) 

-(0.47±0.15)(z-3) 

+ (0.63±0.10)(M 1450 + 26), (10) 



where $(z, < M1450) is the number density of quasars 
brighter than M1450 at redshift z in units of Mpc~ 3 and 
M1450 is the absolute magnitude of the quasar continuum 
in the rest frame at 1450A. Using this equation, we esti- 
mate the number density of quasars, log $(5.7, < —22.0) ~ 
—6.06 ± 0.60 and the number of quasars in our survey vol- 
ume, Aquasar ~ 0.16*q!i2j where M1450 = —22.0 corre- 
sponds to z' w 26.8 and EW w 70 A, if NB816 = 25.0. 
In this estimate, we use a A model with f2 ma ttcr = 0.35, 
Oa = 0.65, and H = 65 km s _1 Mpc -1 , following the 
manner of Fan et al (2001). We note that the survey vol- 
ume for quasars with the broad Lya + N v emission lines 
is larger than the value using above estimate. The survey 
volume for such quasars could increase by a factor of <~ 3. 
We also note that the faint end of the optical luminos- 
ity function of the quasars is expected to be much flatter 
than the bright end which we use in the above estimate 
(see Boyle et al. 2000). If this is the case, the expected 
number of quasars in our survey volume is further lowered. 
In conclusion, the expected number of quasars in our LAE 
candidates is only one at most. 

Note that Malhotra et al. (2003) find no X-ray sources 
associated with 49 LAEs at z = 4.5 or 5.7. It is also noted 
that the results of Steidel et al. (2002) show that only 29 
of 988 LBGs at z ~ 3 have AGN-like emission lines; i.e., 
~ 3 %. These results appear to be consistent with our 
estimate given above. 

3.5. Spatial Distribution of Emitter Candidates 

We show the spatial distributions of the LAE candidates 
at z « 5.7. In Figure 7, we plot the spatial distributions 
for our LAE candidates. The 20 LAE candidates appear to 
be located at west of SDSSp J104433.04-012502.2, while 
the low-z emitter candidates are distributed randomly in 



the sky. However, we do not discuss further detail on this 
issue because spectroscopic confirmation has not yet been 
completed. 

4. DISCUSSION 

4.1. Rest-Frame Lya Equivalent Widths of the LAE 

candidates at z « 5.7 

First, we investigate the rest-frame Lya equivalent 
widths, EWo(Lya), of our 20 LAE candidates (see Ta- 
ble 4). In Figure 8, we show a histogram of EWo(Lya) for 
the 20 sources. Eleven of the 20 LAE candidates are unde- 
tected above 3a level in Iz. For these LAE candidates, we 
estimate lower limits assuming Iz = 26.0 (3a). It is shown 
that the rest-frame equivalent widths range from 25A to 
> 120 A. This range is consistent with that of LAEs at 
z > 5 (e.g., Dey et al. 1998; Hu et al. 2002; Dawson ct al. 
2002). 

4.2. Space Density of the LAE Candidates at z w 5.7 

5.7 



Since we have detected the 20 LAE candidates at z 



objects in the volume of 1.8 x 10 h$ 7 Mpc , we obtain 
the space density of the LAE candidates at z s=s 5.7, 
n(Lya) ~ 1.1 x 10~ 4 ft.Q 7 Mpc~ 3 . Because the above es- 
timate needs correction of the probability of detection, we 
have performed a simulation using the IRAF ARTDATA. 
We assume that emitters have little continuum and are 
point sources. We generated 2,000 model emitters for each 
magnitude interval (Am = 0.5 mag). After smoothing 
the model-emitter images to match to the image size of 
our final image data, these model emitters are put into 
the final image data. Then we try to detect them using 
SExtractor and to select as a LAE with the same proce- 
dure as that used in our paper except for the magnitude 
criterion (NB816 < 25.0). The detection completeness 
of the model galaxies is shown in Figure 9 as a function 
of AB magnitude. From this result, the detection com- 
pleteness is nearly constant at 0.75 for NB816 < 25.0 and 
then decreases for NB816 > 25.0. This analysis shows 
that our detection is fairly well done for NB816 < 25.0 
although we have to correct for the detection complete- 
ness factor when we discuss the star formation rate den- 
sity. Our analysis shows that the detection completeness 
is 0.75 even for brighter sources. The reason for this is 
that a large number of foreground objects affect the de- 
tection procedure in the broad band images (i.e., B and 
i?c) even though our NB816 data are deep enough to de- 
tect emitters with NB816 < 25.0. After correcting for the 
detection completeness, a space density of the LAE can- 
didates is n(Lya) ~ 1.5 x 10~ 4 ft.Q 7 Mpc -3 at z « 5.7. 
On the other hand, Rhoads & Malhotra (2001) obtained 
n(Lya) ~ 8.2 x lO^/i 3 , 7 Mpc" 3 at z « 5.7 for their full 
sample (18 LAEs; see also Rhoads et al. 2003) 9 . This den- 
sity is lower by a factor of 2 than ours. More accurately, we 
need compare the space density of the LAE candidate in 
the same luminosity range since our survey limit of Lya lu- 
minosity is higher than that of Rhoads & Malhotra (2001). 
When their LAE candidates fainter than the faintest Lya 

-|-2 

i 7 ergs s 



luminosity in our sample, 5.0 x 10 42 ft. o7 ergs s 1 , are not 

8 Note that the peak wavelength of the Lya+N v emission of SDSSp J104433.04-012502.2 is ~ 8300 A (see Figure 15). Therefore only a 
part of the Lya flux is probed in our NB81 6 imaging, and thus this quasar is not selected as a LAE in our survey. 
9 n(Lya) ~ 5.9 X 10" 4 ftg 7 Mpc~ 3 for their reduced sample (13 LAEs). 
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taken into consideration, the space density of their LAE 
candidates decreases to n(Lya) ~ 4.5 x 10~ 5 /ig 7 Mpc~ 3 . 
This density is lower by a factor of 3 than ours. This differ- 
ence may be partly due to that their value is not corrected 
for the detection completeness. However, this correction 
seems too small to explain the difference. Therefore, it is 
possible that the difference of the observed fields between 
the two studies (i.e., the cosmic variance). 

4.3. Lya Luminosity of the LAE Candidates 

We investigate the Lya luminosities, L(Lya), of the 
z Rs 5.7 candidates. We assume that all the sources have 
a redshift of 5.70 which corresponds to the case that the 
Lya is shifted to the central wavelength of the NB816 fil- 
ter. In Table 4, we give the Lya luminosities for the 20 
LAE candidates. The derived Lya luminosities range from 
w 5 x I0 42 /i 7 ergs s" 1 to w 1.4 x \Q^h~ 2 7 ergs s" 1 . In 
Figure 10, we show the distribution of Lya luminosities for 
our 20 LAE candidates. We also show the distribution of 
Lya luminosities for the full sample (18 LAEs) of Rhoads 
& Malhotra (2001). Our LAE candidates appear to be 
more luminous than the LAE candidates found by Rhoads 
& Malhotra (2001). The limit of the Lya luminosity for 
our survey, iu m (Lya) w 5.0 x 10 42 /iq 2 ergs s _1 , is larger 
than that of theirs, iu m (Lya) « 2.5 x 10 42 /iQ~ 2 ergs s" 1 . 
This difference explains why our LAE candidates tend to 
be more luminous than those found by Rhoads & Malhotra 
(2001). 

4.4. Star Formation Rates of the LAE Candidates 

We estimate the SFRs for the 20 LAE candidates at 
z « 5.7. Using equation (6), together with a relation 

L(Lya) = 8.7L(Ha) (11) 

from Case B recombination theory (Brocklehurst 1971), 
we can estimate the SFRs using the Lya luminosity; 

SFR(Lya) = 9.1 x 10~ 43 L(Lya) M^r" 1 (12) 

where, L(Lya) is in units of ergs s _1 (see Hu et al. 1998). 
The results are given in last column of Table 4. We note 
that the SFR derived here is reduced by the cosmic trans- 
mission. The SFRs range from 5 to 14 K^Mq yr _1 with 
an average of 8.1 ± 0.3fro jMQ yr" 1 . Although these val- 
ues are typical to the Lyman break galaxies at z ~ 3 - 4 
(e.g., Steidel et al. 1999 and references therein), the num- 
ber density of the strong LAEs like our sources is rather 
small. 

We examine whether or not the SFR derived from the 
Lya luminosity is consistent with that derived from the 
UV continuum luminosity for our sample. The observed z' 
magnitude can be converted to a UV continuum luminos- 
ity at A = 1340 A. Using the following relation (Kennicutt 
1998; see also Madau et al. 1998), 

SFR(UV) = 1.4 x 10- 28 £„ M Q yr" 1 , (13) 

where L v is in units of ergs s^ 1 Hz -1 , we estimate the 
SFR based on the rest-frame UV (A = 1340A) continuum 
luminosity for each object. The results are summarized in 
Table 5. Then we compare the two SFRs, SFR(Lya) and 
5.F.R(UV), for each object in Figure 11. It is shown that 



S'.F.R(UV) appears greater than SFR(Lya) for most of 
our LAE candidates. We also examine global ratio of two 
total SFRs, SFR total (Lya) / SFi? to tai(UV)= 0.77+° ;i 3 , 
where 5'Fi? total (Lya) is the sum of SFR(Lya) of all our 
LAE candidates and SFi?totai(UV) is that of SFR(TTV). 
If only 8 LAE candidates above the 2 a detection in the 
z'-band are taken into consideration, this ratio decreases 
to SFR total (Lya) / SFR total (VV)= 0.47±g$. The lower 
values of S'.F.R(Lya) may be due to the effect of the ab- 
sorption by the dust and scattering by the intergalactic 
medium (sec Hu et al. 2002). 

4.5. Star Formation Rate Density at z w 5.7 

It is interesting to estimate the contribution of the star 
formation in the 20 LAE candidates to the co-moving cos- 
mic star formation rate density (SFRD; e.g., Madau et 
al. 1996). Integrating the SFRs of the 20 LAE can- 
didates, and correcting for the detection completeness, 
we obtain the co-moving SFRD for our sources, psfr ~ 
1.2 x 10- 3 ho. 7 M Q yr" 1 Mpc~ 3 . 

In Figure 12, we compare this SFRD with those of previ- 
ous studies compiled by Kodaira et al. (2003). Note that 
the results of all surveys shown in this figure except for 
the LAE surveys are estimated by converting the Ha or 
UV luminosity densities (obtained by integrating the LFs 
from L = to oo) to the SFRDs by using the equation (6) 
or (13), respectively. As shown in this figure, the SFRD 
derived in this study is much smaller than the previous 
estimates based on the LBG surveys at z = 3 - 5 (Madau 
et al. 1998; Steidel et al. 1999; and Iwata et al. 2003) by 
an order of magnitude. 

On the other hand, our result appears consistent with 
those of the previous LAE surveys for redshift range be- 
tween z ~ 3 and z ~ 6 (Cowie & Hu 1998; Kudritzki et 
al. 2000; Fujita et al. 2003a; Ouchi et al. 2003; Rhoads et 
al. 2003; Kodaira et al. 2003) within a factor of 5. Note 
that all the values shown in Figure 12 are re-estimated 
with the adopted cosmology. Therefore, for readers' con- 
venience, we give them in Table 6. 

The SFRDs estimated from the LAE samples given in 
Table 6 provide only lower limits because of the following 
reasons. First, no correction is made for the extinction by 
dust and/or intergalactic gas. To correct for these effects, 
we need much deeper optical and near-infrared observa- 
tions (e.g. Hu et al. 2002). Second, we do not integrate 
the SFRD of LAE candidates assuming a certain luminos- 
ity function from a lower to a upper limit because there is 
no reliable luminosity function for high-z LAEs. To esti- 
mate the reliable luminosity function for high-z LAEs, we 
need much deeper wide-field surveys to obtain a large sam- 
ple of LAEs. Therefore, future careful investigations will 
be absolutely necessary to estimate a more reliable contri- 
bution of LAEs to the cosmic SFRD at high redshift. 

For future consideration, we try to estimate a Lya lumi- 
nosity function (LF) and then an integrated SFRD using 
the Lya LF. First, we assume that the Lya LF has the 
same functional form as that of Ha LF because it is ex- 
pected that the Lya luminosity is correlated to the Ha 
luminosity as given in equation (11). Here we use the Ha 
LF at z « 0.24 derived by Fujita et al. (2003b) which is 
described as a Schechter function with a = —1.53 ±0.15 
and log0* = —2.62 ± 0.34, where (j>* is in units of Mpc -3 ; 
note that their Ha LF has logL* = 41.95 ±0.25, where L* 



6 



Lymana Emitters at z = 5.7 



is in units of ergs s . This parameter is treated as a free 
parameter in the following analysis. 

Since the small luminosity range covered by our LAE 
candidates, we estimate for the Lya LF satisfying the 
following equation instead of fitting the LF; 



4.6. 
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Follow-up Optical Spectroscopy 

Follow-up optical spectroscopy has been made of two ob- 
jects among the 20 LAE candidates; (1) LAE J1044-0130 
(Ajiki et al. 2002), and (2) LAE J1044-0123 (Taniguchi 
et al. 2003); they are identified as No. 18 and No. 17 in 
Table 2, respectively. These observations were made using 
ESI 10 on the Keck II Telescope and FOCAS on the Subaru 
Telescope. The detailed information is given in the above 
two papers. These two objects are confirmed as an LAE 
at z = 5.687 ± 0.002 and z = 5.655 ± 0.002, respectively. 
In Figure 15, we show their optical spectra obtained with 
ESI on the Keck II Telescope together with that of SDSSp 
J104433.04-012502.2 obtained with FOCAS 11 . 

During the course of our ESI spectroscopy of LAE 
J1044— 0130, we also obtained a spectrum of a galaxy seen 
2.4 arcsec northeast of LAE J1044-0130 (see No. 18 in 
Figure 5). We show this spectrum in Figure 16. Two emis- 
sion lines are clearly seen at A ~ 6710 - 6720 A. These 
lines are identified as [O n]A3727 doublet redshifted to 
z = 0.802 ± 0.002. Most LAEs show a single emission line 
in their optical spectra and thus it is often difficult to iden- 
tify an LAE from other low- z interlopers if the spectral 
resolution is not high enough to check the emission-line 
profile in detail (e.g., Stern et al. 2000). As demonstrated 
in Figure 16, high-resolution spectroscopy (e.g., R > 2000 
- 3000) is very useful to omit low-z [O n] emitters. 

4.7. Optical Counterpart Search for the Lyman Limit 
System at z = 5.72 

Finally, we give comments on a possible counterpart of 
the Lyman Limit System (LLS) at z = 5.72 reported by 
Fan et al. (2000). As shown in Figure 7, the nearest 
LAE to the quasar is LAE J1044-0123. As described in 
detail in Taniguchi et al. (2003), this LAE is located at 
a projected distance of 4.45 Mpc from the quasar if its 
cosmological distance is the same as that of the quasar. 
This LAE is too far from the quasar to be identified as a 
counterpart of the LLS. Further the redshift of this LAE 
(z = 5.655) is different from that of the LLS. Although 
we find a faint galaxy with ms(AB) ss 25 is located at 
1.9 arcsec southwest of the quasar, its broad band color 
properties from B to z' suggest that the galaxy is a fore- 
ground galaxy located at a redshift of z <~ 1.5 - 2.5 (Shioya 
et al. 2002). In conclusion, we cannot find any possible 
counterpart of the LLS. 

5. SUMMARY 

We have presented our optical narrowband (NB816: 
A c = 8150 A and A A = 120 A) observations 
of the field surrounding the high-z quasar SDSSp 
J104433.04-012502.2 obtained with Suprime-Cam on the 
8.2 m Subaru Telescope. The survey area is 720 arcmin 2 . 
Our survey probes a volume of 1.8 x 10 5 /i^ 3 Mpc 3 . 

(1) In our survey, we have found 62 emission- line sources 
whose observed emission-line equivalent widths are greater 
than 180 A. From their optical multicolor properties, we 
have identified 20 LAE candidates at z ~ 5.7 (Az ~ 0.1). 
If they were low-z emitters, they would be very luminous 
star-forming galaxies and thus would be easily detected in 

10 Echclle Spectrograph and Imager available on the Keck II telescope (Sheinis et al. 2000). 

"Faint Object Camera And Spectrograph available on the Subaru Telescope (Kashikawa et al. 2002). 



(14) 

where £(£y lm ) is the Lya luminosity density for the LAE 
candidates with L(Lya) > Ln m and Lu m is the complete- 
ness limit of our survey {L\\ m = 7.0 x 10 42 /i^ 2 ergs s -1 
corresponding to the Lya luminosity of the LAE with 
NB816 = 25.0 whose whole flux in the NB816 pass- 
band is assumed to be attributed to the Lya-emission 
line). Adopting the above Schechter function parameters 
(a = —1.53 and log</>* = —2.62), and using £(£i; m ) = 
oo.oo±/.9 x oo/ 33 (, ,ergs s -1 Mpc -3 for our 15 LAE candi- 
dates with L (Lya) > L lim , we find L» = 4.1±0.4x 10 42 /i^ 2 
ergs s _1 . The Lya LF is shown in Figure 13. Note that 
the detection completeness is taken into account in the 
estimate of £(£r im ). 

Using the same procedure described above, we also 
try to estimate a Lya LF at z « 3.4 using the sur- 
vey data of Cowie & Hu (1998). In this estimate, we 
use £(Ai m ) = A./t oo. i x oo/ 93 (/., ergs s -1 Mpc -3 for 
the 9 objects with L(Lya) > L\ im = 2.2 x 10 42 /i^ 2 ergs 
s _1 . It corresponds to the Lya luminosity of the LAE 
with the limiting narrowband magnitude of their survey 
(HDF LRIS field), toab = 25.0, where we assumed that 
the whole flux in their narrowband (AA = 77 A) imag- 
ing is attributed to the Lya-emission line. We then find 
= 3.8 ± 0.8 x 10 42 /i^ 2 ergs s _1 and the Lya LF is also 
shown in Figure 13 together the result for our LAE can- 
didates at z w 5.7. As shown in Figure 13, there appears 
no evolution between z 5.7 and z « 3.4, although the 
analyses given here are tentative. Note that the detection 
completeness of the survey by Cowie & Hu (1998) is not 
taken into account in the above estimate. 

Then, we re-estimate the SFRDs based on our survey 
and Cowie & Hu (1998) integrating the Lya LFs for the 
two surveys derived above. In order to make it pos- 
sible to compare these LF-corrected results with those 
of LBGs, we take account of the difference in the ex- 
tinction between Lya and UV continuum by adopting 
SFR(Lya)/SFR(UV) = O^t"; 13 (see section 4.4; see 
also Hu et al. 2002). As a result, we obtain a corrected 
SFRD atzfa 5.7, p§f£ - x l0- 2 h . 7 M Q yr -1 

Mpc -3 , where the uncertainty includes that in the esti- 
mate of the Ha LF. We also obtain a corrected SFRD for 
the LAEs at z « 3.4, p§f£ - 2.3±?;g x lQ- 2 h . 7 M Q yr -1 
Mpc -3 . These results are shown in Figure 14. Although 
the lower limits of SFRDs at z « 3.4 and z « 5.7 appears 
significantly smaller than those based on the LBG surveys 
(see Figure 12), the LF-corrected SFRDs are comparable 
to those based on the LBG surveys. If this is real, the 
SFRD appears almost constant between z ~ 1 and z ~ 6. 
Future detailed investigations arc absolutely necessary to 
examine this result. 
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our B and Rq images because their predicted apparent 
magnitudes are fairly brighter than our limiting magni- 
tudes. Note also that two of them have been already con- 
firmed from our follow-up optical spectroscopy as LAEs at 
z = 5.655 (Taniguchi ct al. 2003) and z = 5.687 (Ajiki ct 
al. 2002). 

(2) We find that the rest-frame Lya equivalent widths 
of the 20 LAE candidates range from 25 A to > 120 A, 
being consistent with that of LAEs at z > 5 (e.g., Dey et 
al. 1998; Hu et al. 2002; Dawson et al. 2002). 

(3) After correcting for the detection completeness, we 
obtained a space density of the LAE candidates, n(Lya) ~ 
1.5 x 10~ 4 /ig 7 Mpc~ 3 at z w 5.7. This density is higher by 
a factor of 2 than that of the survey for LAEs at z w 5.7 
by Rhoads & Malhotra (2001). 

(4) We investigate the Lya luminosities for the 20 LAE 
candidates. The derived Lya luminosities range from 
w 5 x 10 42 to 1.4 x 10 43 /i^7 ergs s _1 . Our survey probes 
higher- luminosity sources than that of Rhoads & Malhotra 
(2001) since the equivalent width criterion of our survey 
(.EWobs ;> 180 A) is higher than that of Rhoads & Malho- 
tra (2001; EW ohs > 80 A). 

(5) We obtain SFRs of our 20 LAE candidates rang- 
ing from 5 to 14 ft^M yr" 1 with an average of 8.1 ± 
0.3/iQ yM Q yr -1 . Although these values are typical to 
those of the Lyman break galaxies at z ~ 3 - 4, the number 
density of the strong LAEs like our sources appears rather 
small, since the present survey is not deep enough to de- 
tect faint emission line galaxies. We also compare the two 
SFRs, SFR(Lya) and SFR(VV), for the 20 LAE candi- 
dates. It appears that SFR(UV) is higher than SFR(Lya) 



for most of our LAE candidates. This may be due to 
the effect of the absorption by the dust and scattering by 
the intergalactic medium for Lya emission (see Hu et al. 
2002). 

(6) We have estimated the contribution of the star for- 
mation activity in the 20 LAE candidates to the co-moving 
cosmic SFRD (e.g., Madau et al. 1996). The total star for- 
mation rate of our 20 LAE candidates with correction for 
the detection completeness gives psfr ~ 1.2 x lO~ 3 ft-o.7-^0 
y^ 1 Mpc -3 . This value is much smaller than the previ- 
ous estimates by LBG surveys at z = 3 - 5 (Madau et al. 
1998; Steidel et al. 1999; Iwata et al. 2001) by an order 
of magnitude. However, if we integrate the SFRD adopt- 
ing a Lya LF, the SFRD could be as high as those of the 
previous surveys for LBGs at z ~ 3 - 5. 

We would like to thank both the Subaru and Keck Tele- 
scopes staff for their invaluable help, and Dave Sanders, 
Sylvain Veilleux, S. Okamura, Y. Ohyama, S. Oyabu, N. 
Kashikawa, M. Iye, H. Ando, and H. Karoji for encourage- 
ment during the course of this study. We would also thank 
James Rhoads for providing us useful information on the 
LALA survey, M. Ouchi, M. Yagi, and K. Shimasaku for 
useful discussion on the data reduction of Suprime-Cam 
data, and T. Hayashino for his technical help. We also 
thank an anonymous referee for his/her useful comments 
and suggestions. This work was financially supported in 
part by the Ministry of Education, Culture, Sports, Sci- 
ence, and Technology (Nos. 10044052, and 10304013). 
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Fig. 1. — Response curves (filter, optics, atmosphere transmission, and CCD sensitivity are taken into account) of the bands used in our 
observations. 
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Fig. 2. — Final NB816 image (north is up). Our actual survey area is enclosed by solid line. 
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Fig. 3. — Objects detected to the apparent magnitude limit of NB816 = 26 in the NB81 (>-sclcctcd catalog. The objects selected as 
LAE candidates are shown as filled circles, and other emitters arc shown as open circles. The horizontal solid line corresponds to the color 
of Iz — NB816 = 0.9. Objects above this line have strong emission lines with EW a ^ s fs 180 A or greater. The dashed curves show the 
distribution of 3cr error. 
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Fig. 4. — Diagram of objects detected above 3<r in both Rc~, Ic~ an d z'-band image in the NB816 selected catalog between Rq — Iq and 
Ic — z' (right panel). In the left panel, we show the color-color diagrams expected for five typical type of galaxies(see text). The horizontal 
line shows selection criteria of LAE candidates with 7c < 24.8. 
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Fig. 5. — Broad-band and NB816 images (upper panels) and contour maps (middle panels) of our LAE candidates at z fss 5.7. Each box is 
16" on a side. Each circle is 4" radius. The contours show surface brightness of 25, 26 and 27 mag arcsec -2 . The SED is also shown in lower 
panel for each object. 
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Fig. 6. — B and Rq magnitudes of the model galaxies at z = 0.24, 0.63, and 1.19 shown in three extinction cases (E(B — V) gas = 0, 0.3, 
and 0.6). For readers' convenience, the age of each datapoint is shown slightly shifted from the real value (0.01, 0.1 and 1 Gyr); The B-band 
results arc shifted to the left while the Rc-bsmd ones arc to the right. The model galaxies are assumed to the star forming galaxies with a 
constant SFR estimated from the line fluxes of our LAE candidates. The vertical error bars correspond to the range of SFRs shown in Table 
3. The horizontal line shows our selection criteria of the LAE candidates (i.e., B > 26.6 and Rc > 26.6). 
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Fig. 7. — Celestial positions of the 20 LAE candidates at z 5.7, and 42 lower-2 emitter candidates. The boundary of survey area is also 
shown (solid line). The position of SDSSp J104433. 04— 012502.2 is shown by "x". The two LAEs at 2 ss 5.7 confirmed by our spectroscopy 
are marked. 
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Fig. 8. — Histogram of EWo(Lya) for our 20 LAE candidates. The solid histogram is for seven LAE candidates with the continuum detected 
above 3<r (i.e. Iz< 26.0). The broken histogram is for all our LAE candidates, where EWo of LAE candidates with undetected continuum 
are estimated as lower limits assuming Iz = 26.0 (3cr). 
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Fig. 9. — Detection completeness of LAEs as a function of the apparent NB816 magnitude derived from the simulation. The horizontal 
line shows the detection completeness of 0.75. The vertical line shows the detection limit of our survey (NB816 = 25.0). 
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Fig. 10. — Distribution of Lya luminosity is compared to that derived by Rhoads & Malhotra (2001). We note that Lya luminosities for 
their LAE candidates are estimated from the flux densities shown in Fig. 1 of Rhoads & Malhotra (2001). 
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Fig. 11. — Comparison between SFR(Lya) and 5Fi?(UV) for the 20 LAE candidates. 
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Fig. 12. — Star formation rate density (SFRD) at z ~ 5.7 derived in this study (filled diamond) is shown together with the previous 
investigations. The other data sources are Gallego et al. (1996 - filled triangle), Lilly et al. (1996 - stars), Connolly ct al. (1997 - filled 
squares), Cowie & Hu (1998 - open inverse triangle), Madau et al. (1998 - filled circles), Tresse & Maddox (1998, open circle), Treyer et al. 
(1998 - open triangle), Steidel et al. (1999 - filled pentagons), Kudritzki et al. (2000 - filled inverse triangle), Ouchi et al. (2003 - open 
hexagons), Fujita ct al. (2003a - open diamond), Rhoads et al. (2003 - double triangle), Kodaira et al. (2003 - double circle), and Iwata et 
al. (2003 - open square). 
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dotted lines show the complete limit, L lim , of our survey (7.0 X lO 42 /;,^ 2 ergs s _1 ) and that of Cowie & Hu (1998; 2.2 X lO 42 /^ 2 . ergs s~ 4 ). 



O 
Q- 



L 



10" 1 



10" 



10" 



~i 1 1 1 1 1 r 

This survey ♦ 
Cowie & Hu (1998) v 



1 



t 



6 



Fig. 14. — Corrected SFRD at z f» 5.7 derived in this study (filled diamond) is shown together that of Cowie 
triangle). The shaded area shows the distribution of SFRDs based on the LBG surveys (see Figure 12). 



Hu (1998 - open inverse 



14 



Lymana Emitters at z = 5.7 




S 




6680 6700 6720 6740 

Wavelength (A) 

Fig. 16. — Optical spectrum of the foreground neighbor galaxy of LAE J1044-0130 obtained with ESI on Keck II. The [O n]A3727 doublet 
is clearly seen. 
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Table 1 
Journal of observations 



Band 


Obs. Date (UT) 


Tl nt (sec) a 


m lim (AB) b 


FWHM stal (arcsec) c 


B 


2002 February 17 


1680 


26.6 


1.2 


Rc 


2002 February 15, 16 


4800 


26.2 


1.4 


Ic 


2002 February 15, 16 


3360 


25.9 


1.2 


z' 


2002 February 15, 16 


5160 


25.3 


1.2 


NB816 


2002 February 15-17 


36000 


26.0 


0.9 



a Total integration time. 

b The limiting magnitude (3a) within a 2.8 arcsec aperture. 

c The full width at half maximum of stellar objects in the final image 



Table 2 

Photometric properties of the LAE candidates at z ps 5.7 a 



No. a(J2000) <5(J2000) B a R c a I c a NB816 a z' a Iz a 
h m s ° ' " 



1 


10 


43 


31 


.9 


-01 


14 


46 


>27.0 


>26.6 


(26.2) 


24. 


.6 


(25.5) 


26.0 


2 


10 


43 


38 


.6 


-01 


14 


06 


>27.0 


>26.6 


>26.3 


24. 


.5 


(25.5) 


(26.3) 


3 


10 


43 


43 


.0 


-01 


27 


43 


>27.0 


>26.6 


>26.3 


24. 


.7 


>25.7 


>26.5 


4 


10 


43 


44 


.4 


-01 


15 


59 


>27.0 


>26.6 


25.5 


24 


.1 


24.8 


25.3 


5 


10 


43 


44 


.5 


-01 


24 


31 


>27.0 


>26.6 


(25.9) 


24. 


.7 


(25.5) 


25.8 


6 


10 


43 


45 


.6 


-01 


43 


20 


>27.0 


>26.6 


(26.2) 


24 


.3 


(25.4) 


25.9 


7 


10 


43 


48 


.3 


-01 


23 


20 


>27.0 


>26.6 


25.7 


23 


.9 


(25.5) 


25.6 


8 


10 


43 


53 


.0 


-01 


36 


12 


>27.0 


>26.6 


(26.2) 


25. 


.0 


>25.7 


(26.2) 


9 


10 


43 


55 


.0 


-01 


21 


48 


>27.0 


>26.6 


>26.3 


24. 


.7 


(25.5) 


(26.1) 


10 


10 


43 


55 


.3 


-01 


28 


47 


>27.0 


>26.6 


(26.1) 


24. 


.5 


>25.7 


26.0 


11 


10 


43 


55 


.4 


-01 


13 


39 


>27.0 


>26.6 


>26.3 


24. 


.7 


(25.6) 


>26.5 


12 


10 


43 


55 


.5 


-01 


14 


18 


>27.0 


>26.6 


(26.2) 


25. 


.0 


>25.7 


(26.2) 


13 


10 


43 


59 


.0 


-01 


16 


27 


>27.0 


>26.6 


(26.0) 


24. 


.3 


>25.7 


(26.2) 


14 


10 


44 


08 


.5 


-01 


24 


35 


>27.0 


>26.6 


>26.3 


24. 


.2 


>25.7 


>26.0 


15 


10 


44 


10 


.2 


-01 


16 


18 


>27.0 


>26.6 


>26.3 


24. 


.9 


>25.7 


>26.0 


16 


10 


44 


11 


.3 


-01 


25 


12 


>27.0 


>26.6 


(26.1) 


24. 


.7 


>25.7 


>26.0 


17 b 


10 


44 


26 


.6 


-01 


23 


45 


(26.8) 


>26.6 


(25.9) 


24. 


.7 


>25.7 


25.9 


18 c 


10 


44 


32 


.4 


-01 


30 


34 


>27.0 


>26.6 


(26.1) 


24. 


.4 


>25.7 


26.0 


19 


10 


44 


41 


.4 


-01 


38 


30 


>27.0 


>26.6 


>26.3 


25. 


.0 


>25.7 


>26.0 


20 


10 


44 


52 


.9 


-01 


39 


14 


>27.0 


>26.6 


(26.1) 


24. 


.8 


>25.7 


>26.0 



a AB magnitude in a 2.8 arcsec diameter. The magnitudes between the 2 a and 3 a 
detection levels are put in parentheses. 

b LAE J1044-0123 (Taniguchi ct al. 2003). 

C LAE J1044-0130 (Ajiki et al. 2002). 
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Table 3 

The SFRs assuming the LAE candidates are low-z emitters 





Line 


z c a 


V b 


E(B - V) 


(-^min; -^max) 


(SFR m i n , SFR max ) d 

(H 7 Mq yr j 




Ha 


0.24 


0.39 




0.3 
0.6 


(0.018, 0.054) 
(0.030, 0.089) 
(0.049, 0.146) 


(0.015, 0.043) 
(0.024, 0.70) 
(0.039, 0.12) 


[0 


m]A5007 


0.63 


2.62 




0.3 
0.6 


(0.25, 0.73) 
(0.60, 1.8) 
(1.4, 4.2) 


(0.25, 0.73) 
(0.60, 1.8) 
(1.4, 4.2) 


[0 


n]A3727 


1.19 


7.30 




0.3 
0.6 


(1.2, 3.5) 
(4.5, 13) 
(17, 50) 


(1.1,3.1) 
(4.0, 12) 
(15, 44) 



a The central redshift corresponding to the central wavelength of the NB816 filter. 
b The commoving volume covered by the NB816 filter. 

c Thc minimum and maximum line luminosity, if our LAE candidates are the corresponding line 
emitters. 

d Thc minimum and maximum SFR, if our LAE candidates are the corresponding line emitters. 



Table 4 

LYa luminosity and star formation rate for the LAE candidates at z w 5.7 



No. EW £(Lya) a SFR(Lya) h 

(A) (10 43 /i : 2 7 ergs s" 1 ) (ho 2 7 M Q yr" 1 ) 



1 


63± 6 2 * 


0.8 


7 


2 


>80 


1.0 


9 


3 


>52 


0.8 


7 


4 


46±g 


1.2 


11 


5 


4fl±g 


0.7 


6 


6 




1.0 


9 


7 


100±«» 


1.6 


14 


8 


>36 


0.5 


5 


9 


>57 


0.8 


7 


10 


>72 


0.9 


8 


11 


>51 


0.8 


7 


12 


>37 


0.5 


5 


13 


>110 


1.2 


11 


14 


>118 


1.4 


13 


15 


>43 


0.9 


8 


16 


>58 


0.8 


7 


17 


47+90 

41 '-19 


0.6 


6 


18 


90+ 137 

yu -36 


1.0 


9 


19 


>35 


0.8 


7 


20 


>44 


0.7 


6 


\L): 


=S 0.1 X 10 43 /l^7 


ergs s _1 





h a(SFR) « 1 /t^M Q yr" 1 
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Table 5 

uv continuum luminosity and star formation rate for the lae candidates at z ~ 5.7 



No. 


U{z'Y 

(10-30 crgs s -i cm -2 Hz" 1 ) 


-^1340 


SFR(VV) C 
(h^ 7 M Q yr- 1 ) 


SFR(Lya)/SFR(UV) 




(lO 29 /^ 2 . crgs s- 1 Hz" 1 ) 




1 


2.3 


1.2 


17 


0.44 


2 


2.2 


1.2 


16 


0.55 


3 


<1.9 


<1.0 


<14 


>0.52 


4 


4.5 


2.4 


34 


0.32 


5 


2.2 


1.2 


16 


0.36 


6 


2.6 


1.4 


19 


0.50 


7 


2.3 


1.2 


17 


0.82 


8 


<1.9 


<1.0 


<14 


>0.34 


9 


2.2 


1.2 


17 


0.42 


10 


<1.9 


<1.0 


<14 


>0.56 


11 


2.0 


1.1 


15 


0.48 


12 


<1.9 


<1.0 


<14 


>0.34 


13 


<1.9 


<1.0 


<14 


>0.79 


14 


<1.9 


<1.0 


<14 


>0.91 


15 


<1.9 


<1.0 


<14 


>0.56 


16 


<1.9 


<1.0 


<14 


>0.52 


17 


<1.9 


<1.0 


<14 


>0.42 


18 


<1.9 


<1.0 


<14 


>0.67 


19 


<1.9 


<1.0 


<14 


>0.49 


20 


<1.9 


<1.0 


<14 


>0.43 



*<r(f v ) w 0.9 x 10- 30 crgs s" 1 cm- 2 Hz" 1 . 
b The UV continuum luminosity at A = 1340 A. 

c <t(sfr) » 7 /i .7 M © y r_1 
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Table 6 

The SFRDs derived from previous studies 



No. 




^max) 


PSFR 


Reference 13 








(h . 7 M Q yr- 1 Mpc" 3 ) 




The Ha-emitter surveys 


1 


(0.000, 


0.045) 


1.2 x 10" 2 


1 


2 


(0.0, 


0.3) 


2.5 x 10" 2 


2 


The multi color survcys d 


1 


(0.0, 


0.5) 


1.2 x 10" 2 


3 


2 


(0.20, 


0.50) 


1.1 x 10" 2 


4 


3 


(0.50, 


0.75) 


2.2 x 10" 2 


4 


4 


(0.5, 


1.0) 


4.2 x 10" 2 


5 


5 


(0.75, 


1.00) 


4.2 x 10" 2 


4 


G 


(i.o, 


1.5) 


5.8 x 10" 2 


5 


7 


(1-5, 


2.0) 


4.4 x 10" 2 


5 


8 


(2.0, 


3.5) 


2.9 x 10" 2 


6 


9 


(2.80, 


3.28) 


5.0 x 10" 2 


7 


10 


(3.5, 


4.5) 


1.1 x 10" 2 


6 


11 


(3.87, 


4.39) 


4.0 x 10" 2 


7 


12 


(4.4, 


5.3) 


3.5 x 10" 2 


8 


The LAE surveys 


1 


(3.12, 


3.15) 


6.4 x 10" 3 


9 


2 


(3.41, 


3.47) 


4.7 x 10" 3 


10 


3 


(3.60, 


3.83) 


4.1 x 10" 4 


11 


4 


(4.83, 


4.89) 


6.3 x 10" 3c 


12 


5 


(5.65, 


5.75) 


1.2 x 10" 3 


This paper 


6 


(5.68, 


5.80) 


5 x 10" 4 


13 


7 


(6.51, 


6.62) 


5.2 x 10" 4 


14 



The redshift range of each survey. 

b 1. Gallego et al. (1996), 2. Tresse & Madox (1998), 3. 
Treyer et al. (1998), 4. Lilly et al. (1996), 5. Connolly et al. 
(1997), 6. Madau et al. (1998), 7. Stcidcl et al. (1999), 8. 
Iwata et al. (2003), 9. Kudritzki et al. (2000), 10. Cowie & 
Hu (1998), 11. Fujita et al. (2003a); note that p SFR ~ 1 x 
lO~ 3 ho.7M@ yr" 1 Mpc" 3 when we include their marginal and 
unclassified samples, 12. Ouchi et al. (2003), 13. Rhoads et al. 
(2003) for their full sample (18 LAEs); note that p S FR =^3x 
10" 4 /i .7M Q yr" 1 Mpc" 3 for their reduced sample (13 LAEs), 
& 14. Kodaira et al. (2003). 

c Thcsc SFRs are estimated by integrating the extinction- 
corrected Ha LFs from L = to oo with use of equation (6). 

d These SFRs are estimated by integrating the UV LFs from 
L = to oo with use of equation (13); no correction for dust 
extinction is made. 

°Estimated from the UV continuum of LAE candidates. 
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